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ABSTRACT 

We investigate in detail the mass distribution obtained by means of high resolution 
rotation curves of 25 galaxies of different morphological types. The dark matter con- 
tribution to the circular rotation velocity is well-described by resorting to a dark 
component whose density shows an inner core, i.e. a central constant density region. 
We find a very strong correlation between the core radius size Rc and the stellar ex- 
ponential scale length Ro' Rc — 13 (5^^)^'°^ kpc, and between Rc and the galaxy 
dynamical mass at this distance, Mdyn{Rc)- These relationships would not be ex- 
pected if the core radii were the product of an incorrect decomposition procedure, 
or the biased result of wrong or misunderstood observational data. The very strong 
correlation between the dark and luminous scale lengths found here seems to hold also 
for different Hubble types and opens new scenarios for the nature of the dark matter 
in galaxies. 

Key words: galaxies; kinematics and dynamics - galaxies: spiral - dark matter. 



In spite of a large amount of knowledge gathered over the 
past 20 years on the phenomenon of Dark Matter (DM) in 
Galaxies, it is only recently that the attention has been fo- 
cused on its radial density profile. Recall that, according to 
high-resolution N-body simulations in the A-CDM frame- 
work, the dark halo density, characterized by the virial halo 
mass as free parameter, is described by an inner power-law 
cusp (but see Ricotti 2003): 



PNPw(r-) 



(1) 



where and ps are related parameters (Navarro, Frenk and 
White, 1996b, hereafter NFW). The DM contribution to the 
circular velocity can be written as: 



A{x) 



A{c 

where x = cr / R^i 



(2) 



eter, and A{x 

Vvir ~ GM^ir / Rvir, p. 

Si 
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C)) and M^ir = ^SvirPmeanRlir- ^ 



c = Rvir/Ts is the concentration param- 
ln(l + x) — x/{l + x). We also have that 
101/3 pc„t /(ln(l + c)-c/{l + 
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^ Pcrit is the critical density of the universe, S^ir ~ 337 is the 
virial overdensity and pmean is the mean universal density at the 



A number of reliable mass modellings, obtained for spi- 
rals and low surface brightness galaxies (LSBs) has sup- 
ported the early claim (Moore, 1994) that dark halos around 
disk galaxies have a central density distribution much shal- 
lower than the NFW one (e.g. Gentile et al., 2004 and refer- 
ences therein; Weldrake, de Blok & Walter, 2003; de Blok & 
Bosma, 2002). The first studies on elliptical galaxies seem to 
indicate that also these objects share the same phenomenon 
(Gerhard et al., 2001; Borriello et al., 2003 ). 

This "galaxy by galaxy" comparison between the pre- 
dicted A-CDM density distribution and those actually de- 
tected for the dark halos around disk galaxies has been the 
main goal of several published works. However, although the 
study of the discrepancy between data and theoretical pre- 
dictions remains necessary, we believe that time has come 
that we investigate per se the distribution of DM around 
galaxies, independently of the existing cosmological impli- 
cations. In fact, the hot debate on possible falsification of 
A-CDM on galactic scales has carried the research off from 
the original topic of understanding the DM phenomenon in 
virialized systems. Let us stress that a direct knowledge on 
the presence, nature and interaction with baryons of the 
dark galactic component is still very rough and limited, un- 



galaxy's redshift. The values for the numerical constants are valid 
for a A-CDM cosmogony with f2o=0.3 and Q\=0.7. 
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like the complex and refined scenario that theory and simu- 
lations have put forward. 

In the decomposition of rotation curves (RC's), the cir- 
cular velocity V{r) is best-fit modelled as the sum of a bary- 
onic component that includes a stellar and a gaseous disk, 
and a spherical dark halo: 

V\r)^VEir) + V^l,{r) + V^ir) (3) 

where the labels D, gas and H refers to the correspond- 
ing components. It is common to represent the dark halo 
contribution to the circular velocity by means of the 2- 
parameter function describing the circular velocity of a 
pseudo-isothermal (PI) halo: 

Vlpiir) = ^ttGpoRI (l - ^arctan-^) . (4) 

The central density po and the halo density core radius Rc 
are free parameters to be tuned to fit V^{r). 

Persic, Salucci & Stel (1996, PSS hereafter) introduced, 
for the dark halo velocity contribution, an equally simple 
function, 

VluRcir) = (1 - /3)V;V(1 + a^)/{r^ + a^R%t) (5) 

to fit, once added to the disk contribution, the Univer- 
sal Rotation Curve (URC) of Spirals, i.e. the ensemble 
of synthetic RC's, each one derived from a large number 
of individual RC's belonging to objects with a luminosity 
falling inside a fixed range (see Section 4 of PSS for de- 
tails). The two free parameters are: the velocity core ra- 
dius aRopt, and (1 — /3) Vopt, the halo velocity amplitude 
at Ropt-^ For the aim of this work, in the region of in- 
terest, l/3Ropt < r < Ropt, the PI circular velocity is 
indistinguishable from the halo URC of Eq. 5. In fact, 
the transformation law is obtained from the known quan- 
tities Ropt and Vopt, and by setting the PI free parameters 
Rc and po in Eq. 4 at the values Rc = a Ropt /I -4:5 and 
po = [47rGi?c(l - jg^arctan-^)]-i(l -/3)Kpt. Let us no- 
tice that Vh,pi and Vh,urc agree with each other within 
few percent at most, a discrepancy that is completely neg- 
ligible, in view of the large radial variations of V{r) inside 
galaxies and of the large differences among the halo velocity 
profiles predicted by different mass models. 

Moreover, the velocity profiles of Eqs. 0] and |3 once 
their free parameters have been properly set, can mimic in 
the regions where data are available a number of "theoreti- 
cal" rotation curves, not all necessarily consistent with ob- 
servations."^ 

Therefore, with the obvious caveat that suitable data 
are available, by investigating the properties of the density 
distribution of the dark halos around galaxies by means of 
the PI/URC profile we neither assume nor reject the theoret- 
ical A-CDM scenario. On the contrary, by adopting certain 
values for the free parameters of the chosen halo density 
profile in order to match the observations, we will be able 

^ Ropt is the radius encompassing 83% of the total luminosity of 
the galaxy. In the case of a (stellar) exponential thin disk Ropt is 
3.2 times the disk scale length Rjj- 

^ e.g. the NFW circular velocity with = Ropt, in the ra- 
dial range, 0.3 < r/rs < 1.1 is reproduced, within a negligi- 
ble few percent discrepancy, by Vh URci''",'^ = 0.33) and/or by 
VH,pi{r;Rc = 0.2Ropt). 



to probe a density cusp, or a core or even something in be- 
tween. 

In this work, we first obtain the density profile of DM 
halos around a significant number of galaxies of different 
Hubble types with high-quality kinematic data. The intrinsic 
high quality of the data will allow us to disentangle the 
available kinematics in its dark and luminous contributions 
and to constrain the DM halo parameters allowing us to 
investigate the existence of a link between their structural 
properties and those of the ordinary baryonic matter (i.e. 
stars and gas). 

The sample of data used in our analysis is presented 
and discussed in Sect. |21 while in Sect. |3] we describe the 
mass modelling. In Sect.|l|we will look for a relationship 
between the dark and stellar structural parameters. Finally 
the conclusions will be drawn in Sect.|5] 



2 MASS MODELLING FROM THE ROTATION 
CURVES 

It is necessary, before building a sample of objects with the 
mass model properly derived from the kinematics, to dis- 
cuss the various sources of uncertainties relative to the task 
we want to undertake. In fact, inferring the halo density 
profile from a RC is complicated and it is certain that, in 
spite of applying an effective mass modelling method, we 
get ambiguous or incorrect results when the observational 
errors/biases/uncertainties are larger than a certain thresh- 
old. It is easy to anticipate that a main requirement for a 
proper analysis is that the RC's are of high spatial resolution 
(see below) and that they do not show any significant non- 
circular motion: data must allow us to derive, at different 
radii, a suitable estimate of the circular velocity slope. 

In order to detail this delicate point, let us write the 
halo density p(r) as an explicit function of observational 
quantities, namely the circular velocity, its gradient and the 
disk scale length, and of a derived quantity, the disk mass 
(Fall & Efstathiou, 1980): 

p{r) ~ G" V(r)/r'(l + 2V(r)) - GMoRoptAir/RD) (6) 

where V = dlog V/d\og r, Md is the disk mass and A is a 
given function (for details see the original paper). It is easy 
to show that, since in the inner regions of galaxies the two 
terms in the r.h.s. balance themselves, even mild uncertain- 
ties in V(r) trigger off very large uncertainties in p(r), even 
if Md is known. Small uncertainties in the amplitude and 
in the slope of the circular velocity are an absolute must 
for a proper kinematic density modelling of RC's. Moreover, 
a high spatial resolution is also required in order to have a 
sufficient number of independent measurements to beat ran- 
dom errors and to be able to probe the radial constancy of 
the halo density in the central region, or to follow the NFW 
profile inside the cuspy zone. Notice also that a tight rela- 
tionship between the stellar disk scale length Rd and the 
dark halo core radius such as the one found in this paper 
cannot arise spuriously from uncertainties in the determina- 
tion of Rd, as shown in Appendix A. 

We select our sample of rotation curves of disk systems 
by following the PSS and Borriello and Salucci (2001) pre- 
scriptions. The selected RC must satisfy the following qual- 
ity requirements: 1) data extend (at least) out to the optical 
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Table 1. The sample of galaxies. (1) and (2) name and morpho- 
logical type; (3) stellar disk scale- length Rjj (kpc); (4) and (5) 
core radius Rc (kpc) and central density po (10~^ Mq pc""^); 
(6) references: 1: Borriello & Salucci 2001; 2: Swaters et al. 2003; 
3: Weldrake et al. 2003; 4: Bottema & Verheijen 2002; 5: de Blok 
& Bosma 2002; 6: Bolatto et al. 2002; 7: de Blok et al. 2001; 8: 
Dutton et al. 2003; 9: Simon et al. 2003; 10: Gentile et al., 2004; 
*: J. Simon, private communications 
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radius and to the halo core radius (when detected); 2) it is 
smooth and symmetric; 3) it has small internal r.m.s.; 4j it 
has high spatial resolution and a homogeneous radial cov- 
erage along both arms. The last requirements can be sum- 
marised by setting a maximum uncertainty of 3% in the RC 
amplitude and 0.05 in its logarithmic derivative. 

Let us stress that, when the proper mass modelling is 
applied to the kinematics of galaxies of this suitably selected 
sample, nothing prevents us from detecting a halo core ra- 
dius of zero length or with a value uncorrelated to any other 
galaxy quantity, if this is how Nature has arranged. 

In Table 1 we list the sample of the galaxies selected 
within the above terms that we will use later in the analysis 
and their relevant parameters. The sample is not restricted 
to a single morphological type, but includes dwarfs, spirals, 
LSBs and it considers the most reliable halo density distri- 
butions obtained so far. 

The distances of the galaxies have been taken directly 
in the original papers. They are computed by means of 
local calibrators or assuming Hubble flow after correction 
for Galactic rotation and Virgocentric flow, with Ho=75 
km/sec/Mpc. The uncertainties on the distances can be es- 
timated about 0.30 mag, which implies an uncertainty on 
any galaxy logarithmic linear size of 0.12 dex. 

Finally, let us point out that 1) in literature there are 
about 10 other objects for which the kinematics, although 
it fails our selection requirements, still provides a solid evi- 



dence for the presence of a core radius of which it is however 
unable to constrain the size (e.g. Salucci et al, 2003); 2) we 
failed to find objects of high/average quality kinematics that 
clearly prefer the NFW profile. 

3 MASS MODELS 

Let us detail the mass model, and Freeman disk-|- HI disk 
-I- pseudo-isothermal dark halo used to best fit V{r). 

3.1 Baryonic Disks 

It is usually assumed that light traces the stellar mass via 
a radially constant mass-to-light ratio, so that the stellar 
disk contribution is given by {y = r/Ropt)'- 

vEiy) = 14.7 vE,opt (loKo - IiK,)\,sy (7) 

where Vb.opt = VniRopt)- The circular velocity of the gas 
has also been included and derived from the surface density 
of the HI disk. Gaseous helium and metals have been taken 
into account by a simple rescaling factor of the surface den- 
sity of 1.33. The contribution of the gaseous disk to the total 
rotation speed is found to be marginal. Only in the analysis 
of Borriello & Salucci (2001) and of Bolatto et al. (2002) 
the gas has not been included. However, both the authors 
note that its contribution in the explored region would not 
modify any of their results. 

3.2 Dark Halos 

We choose the pseudo-isothermal halo profile PI, adopted 
to model 20/25 of the sample RC's, as the reference halo 
profile of this work 

where po = Pp/(0) and Rc is the PI core radius. We link 
the core radius of this profile with that of the URC pro- 
file, adopted to model the Universal Rotation Curve, for 
the remaining 5 RC's of our sample (those from Borriello 
& Salucci, 2001) and the sample of ellipticals we consider 
in the conclusions, by recalling that (see previous section) 
the two profiles almost coincide when Rc = ai?opt /I -45 (see 
also, for some examples. Gentile et al. 2004). 

The Rc of NGC 4605 in Bolatto et al. (2002) was fitted 
with in the original paper power-law density p oc r''^'^^ . To 
homogenise it with the analysis on all the other galaxies we 
have fitted their best fit profile in terms of Eq.|H]and derived 
Rc=2.6 kpc. 

We remind the reader that the free parameters of our 
fitting procedure are three, namely: V^{Ropt), po (or /3 for 
the URC RCs) and Rc- Finally, let us point out that, for 
most of the mass models in our Sample a maximun disk is 
not assumed, but the disk is found to be dominant at i?_D , 
as result of a fitting. 



4 RESULTS 

A relationship between "halo core radii" and galaxy lumi- 
nosities was hinted in Persic and Salucci (1988) and in an 
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Figure 1. Core radius Rq as a function of the disk scale length 
R]j. Open and filled circles refer to LSB and HSB galaxies, re- 
spectively. The solid line is the least square fit. 
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Figure 2. Core radius Rc as a function of the dynamical mass 
inside Rc- Open and filled circles refer to LSB and HSB galaxies, 
respectively. The solid line is the least square fit. 



early (unrefereed) analysis of galaxy structure by Kormendy 
(1990). Subsequently, PSS found that the synthetic RC's 
define an Universal Rotation Curve whose dark halo com- 
ponent has a core radius related to disk scale lengths and 
I-band magnitudes: Rc = 3.5Rd 10-(0-2±o i) (Mi-2i.9)/2.5 _ 
This result can be considered as a serious indication that 
the DM distribution relates to the properties of the lumi- 
nous matter. It is, however only indicative in that: 1) the 
RC's coaddition procedure and the relative building of the 
URC was optimised to properly derive the radial behaviour 
of the dark-to-luminous mass ratio rather than the inner 
slope of the dark halo density. 2) The relationship we are 
looking for must be primarily searched for in a sample of 
individual objects, in order to relate the relevant quantities 
in each individual object. 

In this Section we investigate the properties of this "new 
entry" in galactic modelling, the DM core length Rc, i.e. the 
scale length inside which the dark matter density varies very 
mildly, and outside which, it starts to decline as a power law. 
Is Rc only a consequence of the fact that a cuspy distribu- 
tion is a poor representation of actual DM density, or it is a 
real physical quantity related to the dynamical properties of 
the galactic systems? Moreover, could core radii arise from 
observational errors, as recently suggested (Swaters et al., 
2003)? 

The values of core radii for the objects in our sample, 
shown in Table 1, help to answer to these questions, reveal- 
ing a tight correlation between Rc and Rd that we show in 
Fig. 1. This is well fitted by: 

log i?c = (1.05 ±0.11) log + (0.33 ±0.04) (9) 
corr = 0.90, r.m.s. = 0.16 dex 

with corr being the correlation number and r.m.s. the least 
square scatter. The relation clearly exists at very high sta- 
tistical confidence, in addition, LSBs and dwarfs follow the 
spirals' distribution around the mean line with no segrega- 
tion. 

As stressed above, Rd is derived by mere photomet- 



ric measurements, while Rc is obtained from the kinemat- 
ics of the regions dominated by the dark matter. A spuri- 
ous link between these two quantities cannot exist in that 
they are measured in separate ways, and the tight corre- 
lation between Rc and Rd cannot be induced by errors 
in the determination of Rd (see Appendix A). Then, if: 1) 
the cored profile did not describe the actual dark matter 
density in galaxies or 2) the kinematics had serious observa- 
tional errors/biases or 3) the dark/luminous decomposition 
was wrongly made, it is very difficult to imagine how such 
a tight relationship could have emerged. 

We find also a second tight correlation involving the 
core radius and the dynamical galaxy mass inside it, 
Mdyn{Rc) = G-W^{Rc)Rc, that we plot in Fig.H 

logRc = (0.47±0.05) log (^M^^IM^ + (1.15±0.07)(10) 

corr = 0.91, r.m.s. = 0.15 dex. 

This is an additional link between an observed quantity, 
V{r) and the core size Rc and there is no conceivable rea- 
son that would explain such a very tight relation in terms of 
observational or modelling biases, while it is very convinc- 
ing that it simply indicates that the bigger the galaxy, the 
bigger the core radius. Notice that the slope of 0.47, very dif- 
ferent from 1 in Eg. 1101 - that relates two variables each one 
linearly dependent on galactic distance - implies that un- 
certainties on the latter affect very little both relationships 
in Eq.Hand in EalTUl 

Finally, in the case of incorrect mass decomposition of 
RC's, we would expect a particular spurious correlation be- 
tween the core radius and the dark-to-luminous mass frac- 
tion (MF) at Rc, in that both concur and compete to 
shape the RC profile (see Persic and Salucci 1990). In Fig. 
121 we plot Rc as a function of this fraction at r = Rc: 
Mh{Rc)/Md{Rc)- The absence of a tight correlation is ev- 
ident, and demonstrates that the Rc values are not affected 
by gross misfits in the mass decomposition and, more im- 
portantly, that the size of the core radii does not seem to 
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Figure 3. Core radius Re as a function of tlie disk to halo mass 
ratio at Rc- Open and filled circles refer to LSB and HSB galaxies, 
respectively. 



be primarily dependent on the fractional amount of DM in 
this region of galaxies. 

Let us further notice that the above lack of a tight cor- 
relation cannot be related to the existence, at STio, of the 
trend between similar quantities, as implied by the results 
in PSS. By recalling that such relationship is found only 
at fixed multiples of disk scale lengths, we point out that 
the scatter of 0.16 dex in the log Rc vs log Rd relationship, 
though statistically small with respect to the length-scale 
range, is however large enough to imply that Rc corresponds 
to very different galacto-centric distances in different galax- 
ies, i.e. it corresponds to distances varying from I.^Rd to 
?>.1Rd. As a result, the DM fractions computed at iZc's are 
actually computed at a very different radii in different ob- 
jects; and we do not expect to find there the correlations we 
found at each fixed multiple oi Rd. 



5 CONCLUSIONS 

There is a mounting evidence that the kinematics of disk 
galaxies can be explained in terms of the standard disk + 
halo components only if the density distribution of the lat- 
ter decreases with radius very slowly from the center out to 
the radius Rc inside which, instead, the density of CDM 
halos is predicted to fall as a power law with an exponent 
between -1 and -2. Here we have investigated this "obser- 
vational" dark matter scale Rc that is absent in the cur- 
rent theory of galaxy formation. We found that it corre- 
lates tightly with the exponential thin disk scale length Rd 
and with the amount of gravitational mass that it encom- 
passes Mdyn{Rc). The high values for the correlation num- 
ber, the smallness of the scatter around the relations and 
especially the fact that the pairs Rc and Rd, as well as Rc 
and Mdyn{Rc), are measured/derived independently from 
each others, ensure that we are dealing with real and physi- 
cal relationships among physical quantities. 



Three different levels of consequences follow. First, the 
claims according to which core radii arise as a consequence 
of 1) serious observational errors, 2) peculiar or biased kine- 
matics, 3) wrong mass modelling that, until today, have been 
by-passed by considering a number of test-cases, can be now 
ruled out. This, on the basis that human/observational error 
hardly relate to intrinsic properties of galaxies. As an ex- 
ample the suggestion that observers have badly missed the 
galaxy center and doing so artificially created a core, can- 
not explain the additional evidence that this core is found 
to correlate with the galaxy disk scale length within the ex- 
traordinarily small scatter of 0.16 dex. 

The second level realises that the core radius, the quan- 
tity that defines the unexplained feature of the distribution 
of the dark matter, correlates with the main quantity con- 
trolling the distribution of the luminous matter as well as 
the total mass within this radius. This indicates that the DM 
density in galaxies has been shaped by a dark-to-luminous 
matter coupling, a very challenging task if the former is col- 
lisionless. However, let us recall that, at least in dwarf sys- 
tems, according to Navarro et al. (1996a), large mass out- 
flows from galaxies, arising from supernovae explosions, can 
modify an initial NFW profile into a cored one, for which 
Rc oc (Md/ Rd)"'^ . From the results of the mass modelling 
given in PSS, we find Md oc R^, in excellent agreement 
with relationship (9). In any case, our result points toward 
some new physics of which it will be a crucial benchmark. 

The third level is the realization that the core radius vs 
disk scale length relationship is Hubble type free, at least 
within a band of 0.16 dex. For the objects in our sample 
normal/dwarf spirals and LSB lay on the same log Rc vs 
log_RD relationship. Moreover, pioneering analysis on the 
DM distribution in ellipticals (Gerhard et al. 2001) finds 
halo core radii of sizes proportional to the half light radii 
(that correspond to 1.67 Rd): by analysing the mass model 
parameters of 21 Giant Ellipticals (shown in Fig 18 Gerhard 
et al. (2001), and made available by the authors) we find: 



log Rc = (1.1 ±0.2) log Rd + (0.1 ±0.4) 



(11) 



in very good agreement with Eq. (9). We caution that the 
sample of this work is quite limited to fully investigate the 
morphological dependence of DM properties and that in 
Gerhard et al.(2001) the data quality is barely sufficient for 
the intrinsically complicated mass modelling. Nevertheless, 
the emerging picture is truly impressive, especially when we 
consider that core radii of galaxies of different Hubble types 
do not correlate with galaxy luminosity: e.g. a similar value 
for Rc is found in a IO^^Lq^^ spiral, in a lO^°I/0,r LSB 
and in a 10^^1/0, r giant elliptical. The core radius seems to 
uniquely relate to the exponential rate with which the stellar 
surface density decreases with radius. 

The next step of this line of research post-A-CDM cri- 
sis would be to thoroughly investigate this DM distribu- 
tion property in galaxies of different size and Hubble Types 
(Salucci et al., 2004) and the brightest galaxies because of 
the slightly non-linear relationship found between Rc and 
the stellar scale length. 
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6 APPENDIX A 

For sake of completeness, we investigate how much obser- 
vational errors in .Rd could bias the derivation of the halo 
density profile and, in particular, the determination of the 
core radius size Rc ■ We checked this by means of the 4 high- 
resolution RC's of spirals with an exponential thin disk stel- 
lar distribution published in Gentile et al. (2004). First, 
we adopted for Rd the photometric values, we best fit the 
RC's and derived the DM halo core radii Re's. We then let 
the _R_d's to vary by an amount up to ±30%, assuming in 
this way that the photometric values were wrong by this 
amount. Finally, we best- fit again the RC's with these new 
Rd values and derived the corresponding new Rc values. 

As a result, we find that values of Rd even mildly dif- 
ferent from the photometric ones could sometime lead to 
unacceptably high x^'s for any value of the free parame- 
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Figure 4. Changes of values of core radius as a function of the 
disk scale length fractional errors. 



the one obtained for Rnitrue), which is ~ 0.53). This is 
in agreement with the fact that for r < Rd , rotation curves 
are largely dominated by the disk that inprints on it its true 
disk scale-length value (see Salucci et al., 2000; Ratnam & 
Salucci, 2000). 

The study of the fractional error in the estimate of core 
radius as a function of the (assumed) error 

Rc(RD(true)) \ ' 

in the disk scale length "^^^triie)"* shown in Fig. |1] 

We realize that an observational error ixi Rd may trig- 
ger an error ARc in the inferred value of Rc we get by 
fitting the same RC with the wrong value oi Rd'- ARc = 
log{Rc{RD (wrong)) /Rc{Rd (true))). We have 2 possibili- 
ties: 

Large errors (> 10%): ARc can be large. However: i) 
in different objects the same (large) error in Rd induces a 
very different ARc: at a fixed \og{RD{wrong) / Roitrue)) 
the former spans values that are different among themselves 
also by 0.5-1 dex (see Figure 4). This is inconsistent with 
the observed scatter of 0.15 dex of the \ogRc vs logi?_D 
relationship; ii) the slopes of the spurious relationship that 
gets created in each object are completely different from the 
observed one of -1-1 (see Figure 4). 

Small (< 10%) errors: in this most likely case (in our 
sample the typical uncertainty is 5-10%^) ARc turns out to 
be related with disk scale lengths (see Figure 4), but only at 
the level of < 0.1 dex, and therefore the photometric errors 
cannot possibly create the observed range of 1.2 dex shown 
by relationship (9). 

We conclude the investigation by stressing that we can 
consider, for the aims of this paper, the (kinematic) measure 
of Rc and the (photometric) measure of Rd as independent. 



^ Two of these RC's meet the requirements in Sect. 2; the other 
two are, however, of a quality sufBcient for the aim of this Ap- 
pendix: notice below that the effect investigated here leads to 
intrinsic variations of halo core radii much larger than those that 
arise from the observational uncertainties of the RC's. 



^ The uncertainty on Rd was given in the original paper only in 
some cases; in the other cases it was estimated by us. 



